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We analyse the CP asymmetry of the B —* (pKs and B rj'Ks processes in general 
supersjTumetric models. We consider both gluino and chargino exchanges in a model 
independent way by using the mass insertion approximation method. We adopt the QCD 
^ ■ factorization method for evaluating the corresponding hadronic matrix elements. We show 



in 



that chromomagnetic type of operator may play an important role in accounting for the 
deviation of the mixing CP asymmetry between B (pKs and B J/ipKs processes 
observed by Belle and BaBar experiments. We also show that due to the different parity 
' in the final states of these processes, their supersymmetric contributions from the R- 

■ sector have an opposite sign, which naturally explain the large deviation between their 

p • asymmetries. 

• 1 Introduction 



" I 

^ , One of the most important tasks for B factory experiments would be to test the Kobayashi- 

Maskawa (KM) ansatz for the flavor CP violation. The flavor CP violation has been 
studied quite a while, however, it is still one of the least tested aspect in the standard 
model (SM). Although it is unlikely that the SM provides the complete description of 
CP violation in nature (e.g. Baryon asymmetry in the universe), it is also very difficult 
to include any additional sources of CP violation beyond the phase in the CKM mixing 
matrix. Stringent constraints on these phases are usually obtained from the experimental 
bounds on the electric dipole moment (EDM) of the neutron, electron and mercury atom. 
Therefore, it remains a challenge for any new physics beyond the SM to give a new source 
of CP violation that may explain possible deviations from the SM results and also avoid 
overproduction of the EDMs. In supersymmetric theories, it has been emphasised [1] 
that there are attractive scenarios where the EDM problem is solved and genuine SUSY 
CP violating effects are found. 

Recently, BaBar and Belle collaborations [2] announced large deviations from the SM 
expectations in the CP asymmetry of B^ <pKs and branching ratio of B^ r]'K^. 
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These discrepancies have been interpreted as possible consequences of new physics (NP) 
beyond the SM [3-13]. For and 5° decays to a CP eigenstate fcp, the time dependent 
CP asymmetries are usually described by rates af^p{t), 

afcP (t) = ^ fcp) - nB°{t) ^ fcp) ^ ^ ^ ^ 

where Cf^p and S'/^^^ represent the parameters of direct and indirect CP violations re- 
spectively, and AMfi^ is the B° eigenstate mass difference. 

In the SM, the decay process of S — > (f)K is dominated by the top quark interme- 
diated penguin diagram, which do not include any CP violating phase. Therefore, the 
CP asymmetry of B ^ J/ipKs and B — > (pKs in SM are caused only by the phase in 
mixing diagram and we expect Sj/^Ks = S^Ks where Sf^,p represents the mixing 
CP asymmetry. The B rj'Kg process is induced by more diagrams since rj' meson 
contains not only ss state but also uu and dd states with the pseudoscalar mixing angle 
9p. Nevertheless, under an assumption that its tree diagram contribution is very small, 
which is indeed the case, one can expect S^Ks — ^v'Ks ^ well. Thus, the series of new 
experimental data surprised us: 

-^j/^i^s = 0.726 ±0.037, (2) 

which agrees quite well with the SM prediction 0.715lo!o45 [14]. However, results of Belle 
on the corresponding sin 2/3 extracted for B^ — > process has changed dramatically 
[15,16] 

S^Ks = 0.50 ± 0.25l[]:[!I (BaBar), 

= 0.06 ±0.33 ±0.09 (Belle), (3) 

where the first errors are statistical and the second systematic, showing now a better 
agreement than before [17, 18]. However, as we can see from Eq.(3), the relative central 
values are still different. BaBar results [15] are more compatible with SM predictions, 
while Belle measurements [16] still show a deviation from the cc measurements of about 
2a. Moreover, the average S^Ks = 0-34 ± 0.20 is quite different from the previous one 
[19], displaying now 1.7a deviation from Eq.(2). 

Furthermore, the most recent measured CP asymmetry in the B'^ v'^s decay is 
found by BaBar [15] and Belle [16] collaborations as 

Srt'Ks = 0.27 ±0.14 ±0.03 (BaBar) 

= 0.65 ±0.18 ±0.04 (Belle), (4) 

with an average S,^'Kg = 0.41 ±0.11, which shows a 2.5a discrepancy from Eq. (2). For 
the previous results see (BaBar) [20] and (Belle) [18]. 
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It is interesting to note that the new resuhs on s-penguin modes from both experiments 
differ from the value extracted from the cc mode (J/ip), BaBar by 2.7a and Belle by 2.4(T 
[15, 16]. At the same time the experiments agree with each other, and even the central 
values are quite close: 

0.42 ±0.10 BaBar, 0.431°;^^ Belle. 

Supersymmetry (SUSY) is one of the most popular candidates for physics beyond the 
SM. In SUSY models there are new sources of CP violation beside the CKM phase [21]. 
In this review we show our attempts to understand all the above experimental data within 
the Supersymmetric models. 

2 The mass insertion approximation 

As mentioned, the SUSY extension of the SM may provide considerable effects to the 
CP violation observables since it contains new CP violating phases and also new flavour 
structures. Thus, SUSY is a natural candidate to resolve the discrepancy among the 
observed mixing CP asymmetries in i?-meson decays. 

In the following, we will perform a model independent analysis by using the mass 
insertion approximation [22] . We start with the minimal supersymmetric standard model 
(MSSM), where a minimal number of super- fields is introduced and R parity is conserved, 
with the following soft SUSY breaking terms 

VsB = mlM^ + ea,{A^^Yi;HlrA^+A%Yf^^^^^^ 
+ 4rAi/«iXe^. - By^HIHl + H.c.) 

- -(m3^^ + m2W^W^" + miSS) , (5) 

where i, j are family indices, a, b arc SU (2) indices, and Cab is the 2x2 fully antisymmetric 
tensor, with ei2 = 1. Moreover, (f)a denotes all the scalar fields of the theory. Although 
in general the parameters /x, B, and rrii can be complex, two of their phases can be 
rotated away. 

The mass insertion approximation is a technique which is developed to include the 
soft SUSY breaking term without specifying the models in behind. In this approximation, 
one adopts a basis where the couplings of the fermion and sfermion are flavour diagonal, 
leaving all the sources of flavour violation inside the off-diagonal terms of the sfermion 
mass matrix. These terms are denoted by (A^^)*-', where A,B— (L, R) and q — u,d. 
The sfermion propagator is then expanded as 

- . iei - m'l - A^,)- . ^ + (6) 
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where 1 is the unit matrix and rh is the average squark mass. The SUSY contributions 
are parameterised in terms of the dimensionless parameters = (A^^)*-'/rri^. This 

method allows to parametrise, in a model independent way, the main sources of flavor 
violations in SUSY models. 

Including the SUSY contribution, the effective Hamiltonian H^^^ for these processes 
can be expressed via the Operator Product Expansion (OPE) as 

I V ^ p=u,c \ i=3 /J 

+ Q^Q} , (7) 

where \p — Vpiy*g, with Vpb the unitary CKM matrix elements satisfying At + Au + Ac = 0, 
and Ci = Ci{iib) are the Wilson coefficients at low energy scale /if, — mi,. 

As emphasised in [9,11], the leading contribution of both gluino and chargino to AB — 
1 processes come from the chromomagnetic penguin operator Og{Og). The corresponding 
Wilson coefficient is given by 



('51l)23 ( ^M3(a;) + 3M4(a;)) + fc)23^ (^M,{x)+3Msix) 



(8) 



and 

= [(<^ll)32 + mLURsg + mL)s2 + WnLUytR^^ ■ (9) 

Here the functions R^g and Rfg are given by 



Rsg = El^^l|'^W^^^M-«(^^)-^fEW^2Xw.,^P^f fx,), 

K = -Ey^lVi^^W^P^t{x^), (10) 
i 

where x^i = m^/m^., Xi = m^./m^, Xi = m^/m^^, and Xij = rn^Jrn^^. The loop 
functions P^^'^^^\x) and also the functions Mj(a;), i = 1,3,4 can be found in Ref.[12]. 
Finally, U and V are the matrices that diagonalize chargino mass matrix. 

It is now clear that the part proportional to LR mass insertions in Cf^ which is en- 
hanced by a factor rrig/mh would give a dominant contribution. Also the part proportional 
to the LL mass insertion in Cgg is enhanced by m-^/ and could also give significant con- 
tribution. 
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3 Can we explain the experimental data of S(j)Kg in 
SUSY? 



Following the parametrisation of the SM and SUSY amplitudes in Ref.[9], S^Ks can be 
written as 

_ sin 2/3 + 2R^ cos 5i2 sin(g^ + 2p) + Rl sin(2g^ + 2/3) 

l + 2i?^cos5i2Cos^^ + i?J ^ ^ 

where = |A^usy/^sm|^ ^ arg(AsusY/^sM)^ ig g^^^^g pj^^gg_ 

We will discuss in the following whether the SUSY contributions can make S(j,Kg 
negative. For rriq = rrig = 500 GeV and adopting the QCD factorization mechanism to 
evaluate the matrix elements, one obtains 



RT^'^I-, ^ {-0-14 X e-^'-'iSt^hs - 127 x e-''-''(5i^hs} + {L ^ R} , (12) 

while in the case of chargino exchange with gaugino mass M2 = 200 GeV, /i = 300 GeV, 
and TTif^ = 150 GeV, we obtain, for tan /3 = 40 

RT^'X. 1.89 X e-°-°^ (<52J32 - 0.11 x e'^^'^^ (5« J32 

+ 0.43 X e-° '^^ (52^31 - 0.02 x e-°-^^ {61^%,. (13) 

From results in Eqs.(12)-(13), it is clear that the largest SUSY effect is provided by 
the gluino and chargino contributions to the chromomagnetic operator which are pro- 
portional to {Sfj{)-23 and {Sli)32 respectively. However, the 6 — > 57 constraints play a 
crucial role in this case. For the above SUSY configurations, the 6 — > 57 decay set the 
following constraints on gluino and chargino contributions, respectively 1(5^^^)23 1 < 0.016 
and 1(5^^)32! < 0.1. Implementing these bounds in Eqs.(12)-(13), we see that gluino 
contribution can achieve larger value for R^ than chargino one. 

We present our numerical results for the gluino and chargino contributions to CP 
asymmetry S^Ks Fig- 1- We plot the CP asymmetry as function of the phase of {S'[ji)23 
for gluino dominated scienario and A.rg[{6''[j^)s2] for the chargino dominanted model. We 
have scanned over the relavant SUSY parameter space, assuming SM central values as 
in table 1. Namely, the average squark mass m, gluino mass rrig. Moreover we require 
that the SUSY spectra satisfy the present experimental lower mass bounds. In particular, 
rrig > 200 GeV, rh > 300 GeV. In addition, we scan over the real and imaginary part of 
the corresponding mass insertions, by requiring that the branching ratio (BR) of 6 ^ 57 
and the — 5° mixing constraints are satisfied at 95% C.L.. Also we have scanned over 
the full range of the QCD factorization parameters pa,h and 4>a,h, We remind here that 
these parameters are taken into account for the (unknown) infrared contributions in the 
hard scattering and annihilation diagrams respectively. 
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Fi gure 1: S^/iKs ^ ^ function of ai'g[('^L/j)23] (left) and arg[(5j^j^)32] (right) with gluino and 
chargino contributions respectively. 



As can be seen from this figure, the gluino contributions proportional to (5^^)23 have 
chances to drive S<s,Ks toward the region of larger and negative values. While in the 
chargino dominated scenario negative values of cannot be achieved. The reason why 
extensive regions of negative values of are excluded here, is only due to the b — >• 
57 constraints. Indeed, the inclusion of {§11)32 mass insertion can generate large and 
negative values of 5*0, by means of chargino contributions to chromomagnetic operator 
Qsg which are enhanced by terms of order m^±/mb. However, contrary to the gluino 
scenario, chargino contributions to Cgg are not enhanced by colour factors. Therefore, 
large enhancements of the Wilson coefficient Csg, leave unavoidablly to the breaking of 6 — >■ 
57 constraints. As shown in Ref.[ll], by scanning over two mass insertion but requiring 
a common SUSY CP violating phase, a sort of fine tuning to escape 6 — S7 constraints 
is always possible, and few points in the negative regions of 5*^ can be approached. 



4 What happened to the B r]'Ks process? 

Although B —>■ (pKs and B — > rj'Ks are very similar processes, the parity of the final 
states can deviate the result. In _B — >• (j)Ks, the contributions from and Ci to the decay 
amplitude are identically the same (with the same sign), while in B ^ v'^s, they have 
sign difference. This can be simply understood by noticing that 

{(j)KsmB) = {<f)KsmB). (14) 

which is due to the invariance of strong interactions under parity transformations, and 
to the fact that initial and final states have same parity. However, in case oi B ^ v'^s 
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Fi gure 2: Srj'Ks function of <M'g[('^^/j)23] (Isft) and arg[((^^^)32] (right) with gluino and 
chargino contributions respectively. 



transition, where the initial and final states have opposite parity, we have 

{v'KslQmQCDF = -{'n'KsmB)QCDF. 



(15) 



As a result, the sign of the RR and RL in the gluino contributions are different for 
B (pKs and B rj'Ks [23]. Using the same SUSY inputs adopted in Eqs. (13), (12). 
For gluino contributions we have 



QCDF, 



-0.07 X e'^-^\6f 



while for chargino exchanges we obtain 

^ 0.95 {6l,h, 
+ 0.21 (5^^)31 - 0.006 X e 



+ 0.07 xe^o-24(4^)23 + 64(5^^23 (16) 



0.025 X e^*°-^^ (5" 



RL) 32 

-i0.19 / i:u \ 

[<^RL)31- 



(17) 



We show our results for gluino and chargino contributions in Fig. 2, where we have 
just extended the same analysis of -B — > (pKs- Same conventions as in Fig. 1 for B —>■ 
(j)Ks have been adopted here. As we can see from these results, there is a depletion 
of the gluino contribution in 5*^/, precisely for the reasons explained above. Negative 
regions are disfavoured, but a minimum of S*^/ ~ can be achieved. Respect the chargino 
contributions, it is clear that it can imply at most a deviation from SM predictions of 
about ±20 %. 



5 On the branching ratio of B ^ v'^s 

In 1997, CLEO collaboration reported an unexpectedly large branching ratio [24] 

Br^-p-^B'^ K^n') = {89tll ± 9) x 10"^ 



(18) 



which is confirmed by Belle [25] and BABAR [26]: 



BELLE = {79tll ± §) x 10~^ (19) 
BABAR = (76.9 ± 3.5 ± 4.4) x 10~^ (20) 

Considering the theoretical prediction by the naive factorisation approximation 

Br''""' {B K-q') ~ 25 x 10~^ (21) 

the experimental data is about factor of three large, thus, there have been various efforts 
to explain this puzzle. On one hand, new physics contributions have been discussed [27]. 
However, the enhancement by new physics contributions through penguin diagrams ends 
up with large branching ratios for all other penguin dominated processes. Therefore, 
one needs a careful treatment to enhance only B — > -q'K process without changing the 
predictions for the other processes. On the other hand, since this kind of large branching 
ratio is observed only in S — > "q'K process, the ghionium contributors which only exist in 
this process have been a very interesting candidate to solve the puzzle [28] [29] though the 
amount of gluonium in rj' is not precisely known [30]. In this section, we discuss the effect 
of our including SUSY contributions to the branching ratios for B (pK and B — >■ r]'K. 
Inclusion of the SUSY contributions modify the branching ratio as: 

^^SM + SUSY ^ ^^SM X [1 + 2 COS 0svsyR + R^] 

where R = \A^^^^\/\A^^\. As we have shown, to achieve a negative value of S^Ks^ we 
need ^susy — — 7'"/2, which suppresses the leading SUSY contribution. On the other hand, 
the phase for B — > rj'K is different from the one for (f)Ks, as is discussed in the previous 
section. In the following, we will analyze the maximum effect one can obtain from SUSY 
contribution to BR{B — > rj'Ks) with taking into account the experimental limits on the 
BR{B (pKs), S^Ks and S^^^Ks- 

In Fig. 3 we plot the CP asymmetry versus the branching ratio for B (pKs and 
B Tj'Ks- We consider the dominant gluino contribution due to (5^^)23 and scan over 
the other parameters as before. One can see from this figure that, in the region of large 
negative the BR{B — > Kscfy is likely to be close to the SM prediction, namely it is 
of order (2 — 5) x 10~^. Larger values for the BR are also possible but correspond to 
S(t>Ks ^ ~ 0.5. In another word, if we consider the central value of BR{B — > Ks(f)) as 
8 X 10""^, it is predicted that S^Ks hkely to he in the range < Sks4> ^ — 0.5. 

Respect to the correlation between Sr^'Kg and BR(B —>■ Kst]'), it is remarkable that 
with just one mass insertion S^/Rs is likely positive and around 0.5 which is quite com- 
pitablc with the experimental results and in this case with large fi ^ m^, it is no longer 
needed to consider LR and RL mass insertions scmiltancously to suppress 7?^/ as explained 
in the previous section. Furthermore, as can be seen from this figure, for Sri'Xs — 0-5 the 
BR{B rj'Ks) can be large as 80 x 10~^, i.e, it is enhanced by gluino contribution to 
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Figure 3: Correlations of S^Kq versus BR(B 4>Ks) (left) and Sr^iKs versus BR{B rj'Ks) 
(right), for gluino contributions with one single mass insertion (5^^)23- 
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Figure 4: As in Fig. 3, but for gluino contributions with two mass insertions (<^f/j)23 and 

(4l)23. 



more than 6 times the SM value and become compitable with the experimental result 
mentioned above. 

In Fig. 4, we present the correlation between S^Ks BR{B — ^ (pKs) and also the 
correlation between Srj'Xg and BR{B rj'Ks). Here we present the gluino contributions 
with two mass insertions (5^^)23 and (5^/2)32- In this case, we can easily see that the 
gluino contribution can saturate simultaneously both S^Ks and BR{B (pKs) within 
their experimental ranges. However, for B v'^Sy as expected the CP asymmetry 
becomes larger and around the sin 2/9 while its branching ratio is diminished. Now it is 
of order (20 — 40) x 10^^ which is smaller than the experimental measurments. 
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Lr)23 ' (5"ll)32 

2 4 6 8 10 12 14 123456789 10 

BR(B->(^Ks)x10^ BR(B->Ti'Ks)x10^ 

Figure 5: As in Fig. 3, but for gluino and chargino contributions with mass insertions (5f/j)23 
and ((5j^^)32 respectively. 

The combination effects from gluino and cliargino on BR{B (pKs) and BR{B 
rj'Ks) are sliown in Fig. 5. We present tfie CP asymmetry versus the branching ratio 
for each process. We scan on the allowed range of the most relevant mass insertions for 
these two contributions: (5^/^)23 and {S^l)32- We also vary the other parameters as in the 
previous figures. The message of this figure is that with both gluino and chargino we can 
easily accommodate the experimental results for the CP asymmetries and the branching 
ratios of BR{B — > (pKs) and BR{B — > rj'Ks)- It is important to stress that the striengent 
bound on (5^^)32 from the experimental limits on BR{B Xs'y) are relaxed when one 
consider both gluino and chargino contributions, which comes with different sign. Now 
some configuration with large tan/? are allowed and therefore chragino can contribute 
significantly to the CP asymmetries S^Ks and S'^'k^- It is also remarkable that in this 
scienario, the value of the branching ratio BR{B i]'Ks) can be of order 60 x 10"^ which 
is compitable with the central value of the experimental results. 

6 Conclusions 

We studied the supersymmetric contributions to the CP asymmetry of -B ^ (pKs and 
B — i> T]'Ks in a model independent way. We found that the observed large discrepancy 
between Sj/^Ks and S'^^s can be explained within some SUSY models with large ((5Li?)23 or 
('^rl)23 mass insertions. We showed that the SUSY contributions of {5rr)22, and ((5rl)23 
to B ^ (pKs and B — > rj'Ks have different signs. Therefore, the current observation, 
S^iKs < Srj'Ks: favours the (5/jj?,rl)23 dominated models. We also discussed the SUSY 
contributions to the branching ratios. We showed that negative Si^Ks and small SUSY 
effect to Br{B —>■ (pK) can be simultaneously achieved. On the other hand, we showed 



(S''lr)23' (5"ll)32 
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that SUSY contribution itself may not solve the puzzle of the large branching ratio of 
B ^ 7]'Ks. 
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